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Abstract 

Interesting data on dNch/drj in Au-Au collisions (t; — — lntan(S/2)) with the centrality 
cuts have been reported by BRAHMS Collaboration. Using the total multiplicity A^ch ~ 
we find that there are scaling phenomena among (A^ch) ^dNch/drj = dn/drj 
with difi^erent centrality cuts at y'sjvjv = 130 GeV and 200 GeV, respectively. To explain 
these scaling behaviors of dn/drj, we consider the stochastic approach named the Ornstein- 
Uhlenbeck process with two sources. The following Fokker-Planck equation is adopted for the 
present analyses, 



dP{x,t) 



dx 27 dx'^ 



P{x, t) 



where x means the rapidity (y) or pseudo-rapidity (?;). 7 and are the evolution parameter, 
the frictional coefficient and the variance, respectively. Introducing a variable of Zr — ri /rjims 
iVrms = we explain the dn/dzr distributions in the present approach. Moreover, to 

explain the rapidity (y) distributions from ry distributions at 200 GeV, we have derived the 
formula as 

dn j_i dn 
dy drj ' 



where J"^ = y^M{l + sinh^ 2/)/a/1-|- M sinh^ 1/ with M = 1 + (m/pt)^ . Their data of pion 
and all hadrons are fairly well explained by the O-U process. To compare our approach with 
another one, a phenomenological formula by Eskola et al. is also used in calculations of dn/drj. 



1 Introduction 

Recently interesting data on dNch/drj {rj — — lntan(^?/2)) and (0.5(A^part))~^rf^ch/rf'7|j)=o in Au+Au 
collision at y/sNN = 130 GeV and 200 GeV have been reported by BRAHMS Collaboration [1, 2]. 
((A^part) and A^ch niean the numbers of participants (nuclei) and charged particles produced in 
collisions, respectively.) Very recently the BRAHMS Collaboration has reported preliminary data 
on rapidity (y) distribution at 200 GeV in Ref. [3]. We are interested in theoretical analyses of 
these data. 

On the other hand, in Refs. [4, 5] we have investigated the property of 77 scaling of {Nch)^^dNc\i/dr] 
dn/dr] by PHOBOS Collaboration and found that the ij scaling holds. As a possible theoretical 
approach, we have adopted the stochastic theory named the Ornstein-Uhlenbeck (O-U) process 
with two sources at ij/max — ^i^i^/sNN / ttin) ■ In this paper, we would like to analyses data [1, 2, 3] 
by the stochastic approach in terms of the pseudo-rapidity and/or rapidity variables. 

The approach named the O-U process is described by the following Fokker-Planck equation, 
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Figure 1: (a) Initial distribution of Eq. (2). (b) Final distribution at t = finite, (c) Evolution of 
Eq. (2). 



where i, 7 and arc the evolution parameter, the frictional coefhcient and the variance, respec- 
tively ^. Assuming two sources at ij/max — ^i^{y/sNN /ttin) at ^ = and P{y, 0) = 0.5[6{y + 
2/max) + S{y ^ J/max)], wc obtain the following distribution function for dn/dr] (assuming y ^ f]) 
using the probability density P{y, t)[6, 7, 8, 9] 



P{y, ymax, t) 



1 



+ exp 



exp 
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{y-yn 
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where V^{t) = {a"^ /2'y)p with p = 1 — e~'^^*. The physical picture of Eq. (2) with the assumption 
of ?/ « ?7 are shown in Fig. 1. In our approach, it is assumed that A^ch/2 particles are created at 
±J/max at t — 0. Then these A'ch = (A^ch/2 + Nch/2) particles are evolved according to Eq. (2). 
It is worthwhile to mention that a similar approach for the proton spectra has been given in Ref. 
[11]. 

The contents of the present paper are organized as follows. In Sec. 11 i] scaling of BRAHMS 
Collaboration is investigated. In Sec. Ill Analyses of r] distribution by means of Eq. (2) are 
performed. The physical meaning of evolution parameter jt with the frictional coefficient is also 
considered. In Sec. IV Zr = rj/rjrms {Vims = yijf)) scaling is considered. In Sec. V Analysis of y 
distribution derived dn/dr] distribution is presented. In the final section concluding remarks are 
given. 

^The equivalent Langevin stochastic equation with the white noise /w(t) is given as 

dy 



dt 



-7J/ + /w(t). 
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Figure 2: (a) A set of dn/dr] 's with different centrality cuts at ^ysNN = 130 GeV. Each symbols 
have error-bars of about 8 ^ 10% of the magnitude, (b) dn/drj with different centraHty cuts at 
^/sNN = 200 GeV. About the error bars, the situation is the same as (a). 



2 Analysis of rj scaling of dn/drj by BRAHMS Collabora- 
tion 

First of all, we consider the problem on r] scaling in Fig. 2, plotting the data of dn/dr] at 130 GeV 
and 200 GeV. The t] scaling seems to be held. These distributions show 



dn 
dr] 



c (constant). 



(3) 



77=0 



Moreover, we examine the intercept at 77 = 0. Authors of Ref. [11], WA98 Collaboration, noticed 
that the intercepts divided by(0.5(A^part)) should be described by the power-like law, as 



(0.5(iVpart)) 



_1 diVch 



drj 



= ^(A^part)" , 



(4) 



ri=0 



provided that the participants (nuclei) have lost memory and every participant contribute a similar 
amount of energy to particle production in collisions. Actually it can be said that the power-like 
law holds, as is seen in Fig 3. Sec Tables 1 and 2. This physical picture with Eq. (4) indirectly 
supports the availability of the stochastic approach. Combining Eqs. (3) and (4), we have the 
following relations 



1 rfiVch 

iVch dr] 



(5) 



ri=0 



„Sp _ 0-5(iVpart) ^/^j ,Q 

^(^part) , (6) 

where the suffix "Sp" means the scmi-phenomenological formula. Comparisons between Eqs. (5) 
and (6) with A and a in Fig. 3 are shown in Tables 1 and 2. 

As is seen in Tables 1 and 2, the intercept at = is fairly well explained by the semi- 
phenomenological expression, Eq. (6). This implies that the stochastic approach may be available, 
because the participants lost their memory in collision. 



3 Analyses of data by Eq. (2) 

Using the 0-U process with two sources, Eq. (2), we have analyzed the data. The results at 
^/SNN = 130 GeV and 200 GeV are shown in Figs. 4 and 5, and Tables 3 and 4. In our analyses 
we use Eq. (2) the pseudo-rapidity (??) instead of the rapidity (y). As is seen in Tables 3 and 4, 
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Figure 3: (a) Estimation of parameters A and a at y^sTTw = 130 GeV. The method of Unear- 
regression is used. A= 1.64, a = 0.109, and the correlation coefficient (c.c.) is 0.970. A power-like 
law is seen, (b) = 200 GeV. A= 2.09, a = 0.08, and (c.c.)= 0.740. 



Table 1: Empirical examination of Eqs. (5) and (6) at ^/Jnn = 130 GeV. 6ce = 0.013 ^ 0.015 and 
5c. = 0.010 ~ 0.016. 



centrality (%) 


40-50 


30-40 


20-30 


10-20 


5-10 


0-5 


(A^part) 


75 


114 


165 


235 


299 


352 


iVch 


750 ± 60 


1160 ±90 


1720 ± 130 


2470 ± 190 


3180 ±250 


3860 ±430 




0.131 ±5Ce 


0.134±5Ce 


0.138 ±(5ce 


0.141 ±(5Ce 


0.143 ±(5ce 


0.137 ±(5ce 



c^P 0.131 ± 5cs 0.135 ± (5c. 0.137 ±^c. 0.141 ±^c. 0.144 ± Sc. 0.141 ± 5c, 



Table 2: The same as Table 1 but y/s^N = 200 GeV, 5ce = 0.014 ~ 0.016 and 6cs = 0.011 ~ 0.016. 



centrality (%) 


40-50 


30-40 


20-30 


10-20 


5-10 


0-5 


(A^part) 


73 ±8 


114±9 


168 ±9 


239 ± 10 


306 ± 11 


357 ±8 


iVch 


890 ± 70 


1380 ± 110 


2020 ± 160 


2920 ± 230 


3810 ± 300 


4630 ± 370 




0.124 ± (5c, 


0.12G± (5c, 


U.127± (k. 


U.131 ± Sc.,, 


0.133 ± Sc.,, 


U.129± ()Ce 


c"P 


0.121 ±5cs 


0.126 ±Scs 


0.131 ±5cs 


0.133 ±Scs 


0.133 ±5cs 


0.129 ±6cs 
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Table 3: Estimated parameters at y^sjviv = 130 GeV in our analyses by Eq. (2) with two sources. 
Evolution of Eq. (2) is stoppcxl at mininium y^'s. f/max = 4.8. i? = -Avif'^^Z-VcL- 'Auis = \/ (j?^)- 



centrality (%) 


40-50 


20-30 


5-10 


0-5 




789±17 


1775±37 


3273±68 


3952±83 


R 


1.05 


1.03 


1.03 


1.04 


^rms 


2.32 ±0.12 


2.27±0.12 


2.24 ±0.12 


2.21 ±0.12 


p 


0.841± 0.007 


0.858± 0.007 


0.865± 0.007 


0.871± 0.007 


vHt) 


2.79±0.23 


2.80±0.23 


2.64±0.21 


2.56±0.20 


c(Th) 


0.124±0.007 


0.133±0.007 


0.136±0.008 


0.139±0.008 


xVn.d.f. 


0.877/13 


0.434/13 


0.507/13 


0.758/13 




Figure 4: Analyses of dn/dr] at ^/sWn = 130 GeV by Eq. (2). See Table 3. 



R = N^i^ /Nch is always larger than 1 . In the measurements of BRAHMS Collaboration, as the 
observable region is restricted with \r]\ < 4.7, we can conjecture the number of N^'^^^ is always 
3% ~ 7% larger than TVch- 

The different values of in Tables 3 and 4 are attributed to the magnitude of the error bars 
at 130 GeV and 200 GeV. 

The intercepts of dn/dr] at = is explained by the following expression in the 0-U process, 

c(Th) = 1 (exp 

Since our theory is based on the 0-U process, the intercept c'^^'^^ is relating to j/max, the width of 
dn/dr] and the evolution parameter. 

Next we consider physical meaning of the evolution parameter ^t. When we assign the meaning 
of second [sec] to t, j has the dimension of [sec~^]. For the magnitude of the interaction region 
of Au-Au collisions, we assume to be 10 fm. See discussions in Rcf. [12]. See also Tables 5 and 
6. The averaged 7 [fm~^] are almost the same as estimated values from PHOBOS Collaboration 
[13, 14] and ones estimated from the proton spectra at SPS energies in Ref. [10]. 



(±??maxVl - P) 



(7) 
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Table 4: Estimated parameters at ^/snn = 200 GeV in our analyses by Eq. (2) with two sources. 

Evolution of Eq. (2) is stopped at minimum x^'s- Vmax = 5.4. R = N^:'^^^ /Nch, ??rms = ^/ {v'^)^ 
Sp w 0.005 and 6ct = 0.004 ~ 0.005. 



centrality (%) 


40-50 


30-40 


20-30 


10-20 


5-10 


0-5 




955±15 


1477±24 


2158±34 


3101±49 


4034±63 


4881±76 


R 


1.07 


1.07 


1.07 


1.06 


1.06 


1.05 


^rms 


2.41 ±0.08 


2.40 ±0.06 


2.39 ±0.09 


2.37 ±0.08 


2.35 ± 0.08 


2.32 ±0.08 


p 


0.854±(5p 


0.859±(5p 


0M2±5p 


0M6±Sp 


0.871±Sp 


0.878±Sp 




3.169±0.20 


3.17±0.14 


3.15±0.19 


3.16±0.19 


3.10±0.19 


3.08±0.19 


c(Th) 


0.115±(5ct 


0.117±fct 


0.118±(5ct 


0.121±(5ct 


0.123±(5ct 


0.128±(5ct 


xVn.d.f. 


7.2/33 


5.2/33 


4.3/33 


5.4/33 


4.9/33 


5.1/33 
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Figure 5: The same as Fig. 4, but 200 GeV. See Table 4. 



Table 5: Values of 7 and at ^/snn = 130 GcV provided that t w 3.3 x 10 sec. 



centrality (%) 


40 50 


20 30 


5 10 


5 


average 


7 [fm^"'"] 


0.092 


0.098 


0.100 


0.102 


0.098 


[fm-i] 


0.560 


0.601 


0.648 


0.656 


0.616 


aV7 


6.09 


6.15 


6.47 


6.41 


6.28 
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Table 6: Values of 7 and at ^/snn = 200 GeV provided that t rs 3.3 x 10 sec. 



centrality (%) 


40-50 


30-40 


20-30 


10-20 


5 10 


0-5 


average 


7 [fm-1 


0.096 


0.098 


0.100 


0.101 


0.102 


0.105 


0.100 


(7^ [fm-1] 


0.714 


0.723 


0.724 


0.733 


0.729 


0.738 


0.727 


<7V7 


7.42 


7.38 


7.81 


7.30 


7.11 


7.02 


7.26 



The Zr = T]/r]rnis Scaling 



To investigate the Zr = jy/jyrms scaling which has been proposed in Ref. [4], we use ryrms = \/ {ff') = 
if'dn/d'q at ^snn = 130 GeV and 200 GeV. We can consider the following formula with Zr- 



dn dn . , , 

Vrnis~i — ~i ~ J\^r ~ V/Vrms)- 

drj dZr 

The right hand side with multiplying r/rms is obtained from Eq (2), as 



(8) 



dn 

dZr 



1 



cxp 



+ exp 



-7t^2 



2T/2(t) 



{Zr 



(9) 



where z^ax = ??max/(??rms> and V^{t) = (^/ims) IS the averaged quantity in the set of 

data. In concrete analyses of data, V^{t) and p arc treated as the free parameters. The z,. scaling 
at 130 GeV are compared with that of the hemisphere (0 < 77 < 6) at 200 GeV in Fig. 6 (b). It is 
difficult to distinguish them without the labels of incident energies. The behavior of full space is 
given in Fig. 6 (c). This situation is also observed in analyses of data at 130 GeV and 200 GeV by 
PHOBOS Collaboration [13, 14]. 



5 Rapidity {y) distribution derived from ry distribution 

It is well known that one can usually calculate the 77 distribution from the y distribution. In this 
present study, on the contrary, we consider an inverse problem as follows. First we regard Eq. (2) 
as the correct description of the data, because of small values. Using the following formula we 
can obtain the y distribution ^ as 



dn 

_ ^M(l + sinh2 y)dn 
dy ^/l + Msmh^y drj' 

where M = 1 + The right hand side, dn/drj, is given as 



(10) 



dn 
dr] 



1 



+ exp 



exp 



-7*^2 



2F2(t) 



2^2 (f) 



(11) 



2 B-Pz 



■In 



1 + /pj + sinh^ rj + sinh 77 



1 , P + P: 

— In 

2 p-p. 



+ /p1 + sinh^ rj — sinh 77 

dn dn dr] 
dy dr} dy 



■ tanh 



and 



, where 'r){y) = arcsinh(v'M sinhi/) . 



For dn/drj = (j)/E)dn/dy, we have p/E = cosh 77/ \/ 1 + rn^ /p^ + sinh^ r]. Moreover, we have confirmed that 
r°° (dn/dy)dy = 1 and f°° (dn/dr])dr] = 1. 
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Figure 6: Normalized distribution of dn/dzr with Zr = ??/?7rms scaling and estimated parameters 
using Eq. (9). (a) = 130 GeV, p = 1 - ^-27* = 0.889 ± 0.003, V?{f) = 0.527 ± 0.021 

and xVn-d.f. = 5.4/61. (b) and (c) ^/s]^ = 200 GeV, p=l- e'^^* = 0.865 ± 0.002, V^{t) = 
0.559 ± 0.015 and x^/n.d.f. = 32.1/189. (b) is taken from hemisphere data (0 < ?7 < 6) of 
Fig. 6(c). (c) The full space of dn/dZr- The dotted lines represent the magnitude of error-bars in 
the centrality cut 0-5%. 
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where r]{y) = arcsinh(vMsinhy). From Eq. (2) with the averaged parameters p and V^{t), we 
obtain y distributions at 200 GeV for tt meson and all hadrons (tt^, K"^, p and p). They are 
compared with the data in Rcf. [3] in Fig. 7. The small peak is due to the inverse Jacobian factor. 
Indeed the data at 200 GeV show these behaviors at y w 0, even large error bars. To confirm these 
phenomena, measurements in wider region as well as y w are necessary. 

A phenomenological approach proposed in Ref. [15] (which is named as EKRT) is also shown 
in Fig 7. 

dy^ ' {i + e(-y-yo)/d){i + e(v-yo)/d)' ^ ' 

where cn is the normalization factor ^. j/o = 3.3 and d = 0.65 are parameters * given in Ref. [15]. 
Eq. (12) also reproduces the both data in Fig. 7. From Eq. (12) we can calculate dn/dr] (centrality 
cut 20-30%) at 130 GeV and 200 GeV which is presented in Fig. 8. The coincidences between data 
and theory are very well, when yo and d are treated as free parameters. 



6 Concluding Remarks 

1) We have observed that the behaviors of 77 scaling of dn/drj by BRAHMS Collaboration hold 

fairly well among the various centrality cuts at ^/JWn = 130 GeV and 200 GeV. 

2) To explain those scaling behaviors, we have assumed that dn/dr] is governed by the 0-U stochas- 

tic process with two sources at ±?/inax(— In V^ww/tow)- The intercept of dn/drj at = is 
expressed by Eq. (7). See Tables 3 and 4. The constant c's are reflecting the scaling property 

relating to the 0-U process. 

3) From the evolution parameter and the assumed size of the interaction region of Au-|-Au 

collision (10 fm), we have obtained the following value, 7 « 0.1 £m~^, which is almost the 
same value as that estimated in Ref. [11]. 

4) Prom Fig. 6, it can be said that the Zr scaling holds at 130 GeV and 200 GeV. It is difBcult to 

distinguish them, as compared both data without the labels of incident energies. 

5) Using Eq. (11) with rj distributions at 200 GeV, we have calculated the y distributions which 

explain the data of Ref. [3]. The comparison with diS'erent approach given in Ref. [15] is also 
shown. In a future both approaches can be distinguished by the existence of a projection 
(small peak) at y w 0. 

Finally, it can be concluded that the 0-U process is one of possible explanations for the scaling 
property of dn/dr] at ^JsWn = 130 GeV and 200 GeV by BRAHMS Collaboration [1, 2] as well as 
distributions by PHOBOS Collaboration [14]. 
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^We have estimated the normahzation factor c^j as follows 

(1 + e-w/rf)2 

-. r-r-, —, TTi-dy = 6.68. 

(1 + e(-2'-2'o)/'i)(l + e(2'-2'o)/'') 

^Notice that a similar expression with its symmetrization can be seen in Ref. [16]. A different expression based 
on the fractional Fokker-Planck equation for dn/dy is found in Ref. [17]. Both are proposed for analyses of pp (or 
pp) collisions. 
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Figure 7: p and V'^{t) are adopted from Tables 3 and 4. The averaged parameters p = 0.865 (fixed) 
and V'^{t) = 3.138 (fixed) are used, (a) dn/dy of all tt meson, m/pt = 0.4 (fixed), A^ch = 1503±77. 
(b) dn/dy of all hadrons (tt^, K"^, p and p). m/pt = 0.5 (fixed), A^ch = 3915 ± 234. The dashed 
lines are obtained from Eq. (12) with 1/cn = 0.149. 
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Figure 8: Using Eqs. (11) and (12), we calculate dNch/dr] (ccntrality cut 20-30%) at 130 GeV and 
200 GeV. (a) Dashed line is obtained by yo = 3.3, d = 0.65, N^'^^'^ = 1731 ± 34, m/pt = 0.5 and 
xVn.d.f. = 12.1/15. Solid line is obtained by 0-U process [xVn.d.f. = 0.43/13 from Table 3], (b) 
Dashed line is obtained by yo = 3.3, d = 0.65, Nj^ = 2054±31, m/pt = 0.5 and xVn.d.f. = 134/35. 
Solid line is obtained by 0-U process [x^/n.d.f. = 4.3/33 from Table 4]. When yo and d are 
treated as free parameters, the following sets of parameters are obtained, (c) yo = 3.32, d = 0.83, 
xVn-d.f. = 0.40/13. (d) yo = 3.72, d = 0.83, xV^-d.f. = 4.4/33. 
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